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£N| . In homogeneous cosmologies, quantum geometry effects lead to a resolution of the classical 

singularity without having to invoke special boundary conditions at the singularity or introduce 
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ad-hoc elements such as unphysical matter. The same effects are shown to lead to a resolution 
of the Schwarzschild singularity. The resulting quantum extension of space-time is likely to have 
significant implications to the black hole evaporation process. Similarities and differences with the 
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^ ' situation in quantum geometro dynamics are pointed out. 
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I. INTRODUCTION 



General relativity provides a subtle and powerful interplay between gravity and geometry, 
thereby opening numerous possibilities for novel phenomena. Among the most spectacular 
of the resulting conceptual advances are the predictions that the universe began with a big 
bang and massive stars can end their lives as black holes. In both cases, one encounters 
singularities. Space-time of general relativity literally ends and classical physics comes to a 
halt. 

However, general relativity is incomplete because it ignores quantum effects. It is widely 
believed that quantum gravity effects become significant in the high curvature regions that 
develop before singularities are formed. These are likely to significantly change the space- 
time structure, making the predictions of general relativity unreliable. Hence, real physics 
need not stop at the big bang and black hole singularities. While the classical space-time 
does end there, quantum space-time may well extend beyond. 

Attempts at extending physics beyond the big-bang singularity date back at least to 
the seventies when mini-superspaces were introduced and quantum cosmology was born 
(see, e.g., jlj]). These investigations sparked off new developments in different directions 
ranging from foundations of quantum mechanics to the development of WKB methods to 
test semi-classicality in quantum cosmology, to the introduction of novel Euclidean methods 
to calculate the appropriate wave function of the universe. This deepened our understanding 
of issues related to the quantum physics of the universe as a whole However, in these 
approaches the singularity was not generically resolved. A quantum extension of space- 
time either required the introduction of new principles jj, or ad- hoc assumptions, such as 
existence of matter violating energy conditions already at the classical level, or of external 
clocks that remain insensitive to the infinite curvature encountered at singularities. 

Emergence of quantum Riemannian geometry and the associated mathematical tech- 
niques 0, 0, |f| have provided a new approach to revisit the issue of quantum extensions. 
Within the mini-superspaces used in quantum cosmology, the big-bang singularity could 
be resolved without having to introduce external boundary conditions at the big-bang, or 
matter/clocks with unphysical properties 

BBQ. 

In the homogeneous, isotropic mini- 
superspace coupled to a massless scalar field, in particular, an exhaustive analysis can be 
carried out [10]. At the analytic level, one can introduce an appropriate inner product on 
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physical states, define Dirac observables and, using the two, construct semi-classical states. 
One can then use numerical methods to examine the nature of the quantum space-time. If 
one begins with a semi-classical state at large late times ('now') and evolves it back in time, 
it remains semi-classical till one encounters the deep Planck regime near the classical singu- 
larity. In this regime the quantum geometry effects dominate. However, the state becomes 
semi-classical again on the other side; the deep Planck region serves as a quantum bridge be- 
tween two large, classical space-times. The model is too simple to be applied reliably to the 
actual universe we live in. However at a conceptual and mathematical level it demonstrates 
the new possibilities, systematically made available by quantum geometry. 

It is then natural to ask if a similar resolution of singularities occurs also in the context 
of black holes. Now, there is an extensive literature on the nature of black hole singularities 
in the classical theory. In particular, detailed mathematical analysis of spherical collapse 
of uncharged matter was carried out by Christodoulou, Dafermos and others. It has led to 
the expectation that, for black holes formed through gravitational collapse, the singularity 



would be generically space-like (see, e.g., JJ| and references therein). Therefore, it is natural 
to focus on this case first and ask for the nature of quantum space-time that would result 
via their resolution. 1 If the singularity is resolved and if the quantum geometry in the deep 
Planck regime again serves as a bridge to a large classical region beyond, there would be no 
information loss in the black hole formation and evaporation process (see also Q, Q|)- 
Pure states in the distant past would evolve to pure states in the distant future and one 
would have a space-time description of the entire process in a quantum mechanical setting. 
Thus, a detailed analysis of the fate of black hole singularities is of considerable importance 
also for the fundamental issue of whether the standard unitary evolution of quantum physics 
has to be modified in the setting of black holes. 

The purpose of this note is to initiate this investigation using quantum geometry, in 
the setting of connection dynamics. We will focus just on the mini-superspace that is 
appropriate for describing the geometry interior to the horizon of a Schwarzschild black 



1 Indeed, one does not expect quantum gravity effects to resolve all singularities. If, for example, they led 
to a resolution of the time-like singularity of the negative energy Schwarzschild space-time, energy would 
be unbounded below in quantum gravity and the theory would have unphysical features [l2|. Rather, one 
expects that there would be no physical states in quantum gravity that would resemble negative energy 
Schwarschild geometry at large distances, whence the issue of 'resolution' would simply not arise. 
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hole. Since this mini-superspace is spatially homogeneous — of Kantowski-Sachs type- 
one can take over the techniques that have been developed in the setting of homogeneous 
cosmologies |S]. Our main result is that the quantum scalar constraint is indeed such that 
the singularity is resolved. The salient features of our analysis are as follows. First, we 
use a self-adjoint Hamiltonian constraint. Therefore, our analysis can serve as the point of 
departure for the construction of the physical Hilbert space either via the group averaging 
procedure [kJ or via deparametrization of the theory as in Second, we spell out the 
symmetry reduction procedure that leads to the Kantowski-Sachs mini-superspace within 
connection dynamics. We will see in sections IIHI and HVl that the structure of the resulting 
phase space is important to the issue of singularity resolution. 2 In particular, this space is 
an extension of the phase space used in geometrodynamics in that one allows for degenerate 
triads (and hence 3-metrics). Thanks to this enlargement, points of the reduced phase space 
corresponding to the singularity do not constitute a boundary, whence the support of the 
wave function 'beyond singularity' can be interpreted geometrically. However, we emphasize 
that the results of this note constitute only initial steps for a more complete theory, e.g., 
along the lines of jlfl| . 

Our discussion is organized as follows. In section [H] we discuss the classical theory of the 
Kantowski-Sachs mini-superspace using connection-dynamics. In section ITTTl we analyze the 
kinematics of its quantum theory and in section HVl we present quantum dynamics. Section 
13 provides a brief summary and directions for future work. 

II. CLASSICAL THEORY 

This section is divided into two parts. In the first, we carry out the Kantowski-Sachs 
symmetry reduction of vacuum general relativity and in the second we discuss the structure 
of the resulting phase space. 



2 A q uantization of the Kantowski-Sachs minisuperspace is available in the geometrodynamical framework, 



A q 





Unfortunately, it does not shed light on singularity resolution because the location of the singularity 
in the minisuperspace was misidentified. In the notation used in 17], the singularity lies at (a = oo, b = 0) 
— not at (a = 0, b = 0) as assumed there. Therefore, neither is the inverse volume operator (17) of [17J 
bounded or well-defined at the singularity nor does the discrete equation (23) of enable one to evolve 
across the singularity. 
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A. Symmetry reduction 

The portion of the space-time interior to the horizon of a Schwarzschild black hole can 
be naturally foliated by 3-manifolds M with topology Ix § 2 such that the geometry on 
each slice is invariant under the Kantowski-Sachs symmetry group G = R x SO (3). In 
the discussion of this model we closely follow |9j where the mathematical structure of the 
simpler, isotropic model is studied in detail. 

As in any homogeneous situation, M acquires an equivalence class of positive definite 
metrics related by an overall constant. We fix one such metric °q a \, as well as an orthonormal 
triad °e° and co-triad %j % al compatible with it. With the given symmetry, we choose the 
metric on the orbits of the SO(3) action to be the unit 2-sphere metric, d$ 2 + sin 2 fidip 2 
in polar coordinates $, (p. Locally, the corresponding co-triad elements d$, and siniMy? 
can be completed to an orthonormal co-triad by a third form a. Thanks to the Maurer- 
Cartan relations for the symmetry group, a must be closed and therefore exact on the given 
topology, whence it defines a third adapted coordinate x via a = dx. In these coordinates 
the background metric has line element 

ds 2 D = dx 2 + d$ 2 + sin 2 $d<p 2 

and determinant °q = sin 2 

In connection dynamics, the canonically conjugate pair consists of fields (A l a , E?) of an 
SU(2) connection 1-form A l a and a (possibly degenerate) triad Ef of density weight 1 on M. 
This pair is invariant under the symmetry group if it satisfies: 

C^A = DA and C^E = [E, A] (1) 

for any symmetry vector field £ a and some generator A* of local, SU(2) gauge transformations 
(which may depend on £ a ). One can verify that any such symmetric pair is gauge equivalent 
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to fields of the form: 3 



A = cr 3 dx + (art + br^d'd + (— br\ + ar 2 ) sin^dy? + r 3 cos fidip (2) 

d d 

E = p c r 3 sin^^: + (pan + PbT 2 ) sin ^^ + (-Ptn + p^)q^ ■ ( 3 ) 

where a, 6, c are real constants and r J is the standard basis in su(2), satisfying [r», Tj] = Tfc. 
However, this does not exhaust the freedom to perform local SU(2) gauge transformations 
entirely: There is still the freedom to perform a global SU(2) transformation along r 3 which 
rotates the 'vectors' (a, 6) and (p a ,Pb)- We will return to this freedom in the next sub-section. 
From the invariant triad density (J3J we can derive the corresponding co-triad uj: 



uj = uj c r 3 dx + {uJaTi + uJ b T 2 )d'& + {-uJbTi + ^ a T 2 ) sin^dy? (4) 



where 



sgnp c 0>a + pi y/\pT\Pb , ViKlPa , K v 

V\Pc\ VPa+Pb VPa+Pb 

The co-triad in turn determines its spin-connection, the su(2)-valued 1-form T. As one might 
expect from the structure of the triad, T turns out to be the standard 'magnetic monopole' 
connection: 

T = cos^t^ dip . (6) 
Consequently, the su(2)-valued extrinsic curvature 1-form K is given by: 

■yK : = A — T = cr 3 dx + (art + br 2 )d'd + (—bri + ar 2 ) sin$dy? . (7) 

Finally the curvature F a b of A a can be easily computed. Its dual, the su(2)-valued vector 

density B a := \rf hc F bc has the same invariant form as E a : 

d d d 

B = (a 2 + b 2 - l)r 3 sin$— + (an + br 2 )c sin + (-br x + ar 2 )c— (8) 

ox ov dip 



These expressions will be used in sections IIHI and |TV] 



3 More precisely, invariant connections carry a non-negative, integer-valued topological charge which 

we have set equal to one in the above expressions. (In general, it would multiply the last contribution to 
see, e.g., |2jJ.) Each value of this charge gives rise to an independent sector of invariant connections. How- 
ever, only the sector used here allows non-degenerate triads invariant under Q with the same local gauge 
transformation in A and E. That A in J5J satisfies (Q) can be verified as follows. If we denote the three 
generators of the SO (3) action on M by X = smLpd$ + cot i? cos ipd v , Y = — cos ipd^ + cot 1? sin </3cL, and 
Z = d v , and the generator of the R-action by t a , thenCxA — D(cosf/ sinz?)T 3 , LyA — D(siny>/sini?)T3, 
CzA = and CtA = 0. Invariance of E in follows analogously, or by noting that Ef8A l a defines a 
gauge invariant and (spatially) constant 1-form on phase space. 
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B. The reduced phase space 



The symmetry reduction procedure led us to a 6-dimensional reduced phase space T with 
coordinates (a, b, c; p a ,Pb,Pc)- Let us begin by computing the symplectic structure fl on it. 
The basic idea, as in all symmetry reductions, is to pull-back the symplectic structure of 
the full theory to the symmetry reduced phase space. However, since the full symplectic 
structure involves an integral over M and since the fields of interest are homogeneous also 
in the non-compact M-direction, as usual, we are led to consider only a finite interval X in 
the R direction. Let the length of this interval (w.r.t. the fiducial metric °q a b) be L . Then, 
the symplectic structure on the reduced phase space is given by: 

Cl = (2da A dp a + 2db A dp b + dc A dp c ) (9) 
27G 

where 7 is the Barbero-Immirzi parameter and G the gravitational constant. (The volume 
of the § 2 orbits of the SO (3) action, defined by the fiducial metric, is the standard one, An.) 

Because of the form (j2J) and (jSJ) of invariant connections and triads, the vector constraint 
is automatically satisfied. However, because of the residual, global SU(2) gauge freedom 
mentioned in section Hi Al the Gauss constraint is not. Inserting the invariant connection 
and triad into the Gauss constraint we obtain 

C G auss = ap b -bp a = (10) 

which generates simultaneous rotations of the pairs (a, b) and (p a ,Pb)- Thus, only the 'norms' 
\/ a 2 + b 2 , \fp 2 + p\ and the 'scalar product' (ap a + bp b ) are gauge invariant. We will fix this 
gauge freedom as follows. If (a, b) = 0, we rotate the triad components such that p a = 0. 
Otherwise we rotate (a, b) such that a = which implies b 7^ 0. Then, Ccauss = implies 
p a = 0. There is still a residual gauge freedom U b : (b,p b ) — > (—6, — p b ) which changes the 
signs of b and p b simultaneously. This is just the parity transformation in the p b variable. 
One can either retain this freedom and ensure that all the final constructions are invariant 
with respect to this '6-reflection' or eliminate it by a gauge choice such as p b > 0. This 
gauge choice turns out to be particularly convenient for classical dynamics. In the quantum 
theory, on the other hand, it is more natural to retain the freedom at first and then ask that 
physical states be invariant under the parity operator fl b implementing this transformation. 
Therefore, we will allow both possibilities. In either case the 4-dimensional phase space 
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carries a global chart (b,c;pb,p c ) (where pb > if the '6-parity gauge' is fixed.) These 
coordinates are subject to the scalar or the Hamiltonian constraint, discussed below. 

We first note that these variables and, because of its explicit dependence on L a , the sym- 
plectic structure Cl depend on the fiducial metric. It is convenient to remove this dependence 
by rescaling the variables in a manner that is motivated by their scaling properties, 

(b, c) := (6, L c) , (p b ,p c ) := (L p b ,p c ) . (11) 

We now have a one-to-one parametrization of the gauge fixed phase space T by (b,p b , c, p c ) 
with symplectic structure 

tl = -^—(2dbAdp b + dcAdp c ). (12) 
Z7G- 

For later purposes let us express the volume of our elementary cell Ix§ 2 and areas bounded 
by our preferred family of curves as functions on the reduced phase space. The volume is 
given by: 

V = J d 3 xy/\detE\ = ATxL ^\p\\\p h \ = ^^[p^Wp^ . (13) 

The three surfaces S X}(p , S Xj # and S& tV of interest are respectively bounded by the interval X 
and the equator, X and a great circle along a longitude, and the equator and a longitude (so 
that Sd iV> forms a quarter of the sphere S 2 ). Their areas are given by: 

A x ^ = 2n\p b \, A x>(p = 2n\p b \ : and = n\p c \ . (14) 

Finally, let us consider the Hamiltonian constraint and classical dynamics. On the reduced 
phase space T, the constraint functional of the full theory 

C Ham = J dx'Ne- 1 [e ijk Fi h E«E b k - 2(1 + ^K^K^E^ , (15) 



where e := y^jdet E \ sgn(det E), reduces to 



87rA^ Sgn Pc r 2 2 - 2 n , 

C H am = 2 j=^- [{b +7 )Pb + 2cp c bpb\ , (16) 

7 V\Pc\Pb 



where, as is usual in the homogeneous models, we have chosen a constant lapse function 

WnGb 



N. To simplify the equations of motion, it is convenient 4 to choose iV = Zfg^l_^jfd m Then 



4 With this choice, the points b = or p c = have to be excised in the discussion of dynamics. The explicit 
form of the solutions below shows that this is not overly restrictive since on classical solutions these points 
correspond to the singularity or the horizon. 
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the Hamiltonian constraint becomes Cuam '■= — (27G) l [{b 2 + 7 2 )pb/& + 2cp c ], yielding the 
following equations of motion: 

b = -\{b + 1 2 b- 1 ), p h = \{pi > -l 2 b- 2 p b ) 

c = -2c, p c = 2p c , 

where the 'dot' denotes derivative with respect to the affine parameter of the Hamiltonian 
vector field. The equations for (c,p c ) can easily be solved by c(T) = c e~ 2T and p c (T) = 

,(0) 2T 



p c e 2T '. For the other components we obtain b(T) = ±7\/ e~( T ~ To ) — 1 for 6 and, using this 
result, Pb{T) = pf* 'v 'e T+T ° — e 2T . If we introduce the new time parameter t := e T and the 
constant m := |e T °, this gives 

b{t) = ± 7v /(2m-t)/t, p 6 (t) = y/t(2m - 1) (17) 

c(t) = T7^6 *~ 2 , Pc(*) = ± * 2 (18) 

where we have fixed the multiplicative constant in the expressions of p c {t) and c(t), respec- 
tively, by requiring that \p c \ be the geometric radius of the 2-sphere orbits of the SO(3) 
symmetry and by using the Hamiltonian constraint. Projections to the p^ — p c plane of 
typical trajectories are shown in figure H It is clear from the figure that p c can be taken to 
be the 'internal time parameter'. In section HV Al this interpretation will let us regard the 
quantum Hamiltonian constraint as an 'evolution equation' with respect to an internal time 
parameter defined by the eigenvalues of the operator p c . 

Note that we have a 2-parameter family of solutions, labeled by m and pf\ as one would 
expect from the fact that the reduced phase space is 2-dimensional. However, in a space- 
time description we expect only a 1-parameter family, labeled by m. This discrepancy can 
be traced back to the standard tension between the Hamiltonian and space-time notions of 



gauge (for a detailed discussion in the context of Bianchi models, see 21], and in spheri- 
cally symmetric models, 22, 13]). In the Hamiltonian description, rescalings of p^ do not 
correspond to gauge because they are not generated by any constraint. In the space-time 
description, on the other hand, they can be absorbed by a rescaling of the coordinate x. 
Therefore, to make contact with the space-time description, let us fix = 1. Then, the 
above solution to the evolution equations defines a 1-parameter family of 3-metrics: 

/2m \ 

qabif) = (--!) V a xV b x + t 2 (V a $V b $ + sin 2 #V a <pV b <p) (19) 
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FIG. 1: Dynamical trajectories in the pb — p c plane. Each trajectory reaches the maximum value m of pb 
and meets the p c — axis only at the point where pb also vanishes (i.e., at the 'origin'). Solutions related 
by the reflection-symmetry p c — > —p c define the same metric but carry triads of opposite orientation. For 
simplicity the '6-parity gauge' is fixed by requiring pb > 0. 

defining precisely the Schwarzschild solution of mass m, interior to the horizon. It will be 
useful to note that, in this solution p c > 0, p b = at the horizon (t = 2m) and p c = and 
b, c diverge at the singularity. 

We will see in the next section that it is convenient to regard the space Ai spanned by 
p b and p c as the mini-superspace on which quantum wave functions are defined. (In that 
discussion, we will not fix the '6-parity gauge'; so p c and p b take values on the entire real line.) 
Let us therefore interpret various regions in Ai. The triad is non-degenerate everywhere 
except when p b = or p c = 0. However, the two degeneracies are of very different geometric 
and physical origin. Geometrically (i.e., independent of classical dynamics) p c = separates 
two regions with triads of opposite orientations, given by sgndet-E = sgn(p c pl). On the 
other hand, as (@J) (together with our gauge choice which ensures a = p a = 0) shows, the 
orientation would not have changed even if we had extended Ai across p b = 0, allowing 
for negative values for p b . More importantly, while the co-triad (@J) remains smooth and 
only becomes degenerate along p b = 0, it diverges along p c = 0, signaling the classical 
singularity. This suggests that the line p b = corresponds to the horizon and p c = to 
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the singularity. This is confirmed from physical considerations involving classical dynamics. 
Classical solutions can be smoothly extended to the line pb = which is reached for t = 2m. 
The curvature B remains well-defined there (see (jSJ))- The other line, p c = 0, is very 
different. Along classical solutions, the line p c = is approached as t — ► 0, whence the 
components (and invariants) of curvature B also diverge. This line corresponds to the 
singularity and cannot be crossed by any classical dynamical trajectory. Finally, note that 
since the line p c = separates the two regions of Ai, where only the triad orientation is 
reversed, the two get identified in geometrodynamics. Therefore the singularity would lie at 
a boundary of the mini-superspace in geometrodynamics, making the geometrical meaning 
of 'evolution across the singularity' obscure. Furthermore, since the Arnowitt-Deser-Misner 
(ADM) variables are based on covariant metrics (rather than contravariant triad densities), 
as in other homogeneous models one of the ADM variables becomes infinite at the 
classical singularity. This further complicates a discussion of the singularity structure and 
its resolution. 



III. QUANTUM KINEMATICS 

This section is divided into two parts. In the first we introduce the basic quantization 
procedure and in the second we discuss quantum geometry with emphasis on the classical 
singularity. 



A. Basics 

As in the isotropic model , we will follow the general procedure used in the full theory 
0, E|. Thus the elementary configuration variables will be given by holonomies along 
curves in M and the momenta by fluxes of triads along 2-surfaces in M. However, because 
of symmetry reduction, one need not consider all (piecewise analytic) curves or surfaces; 
a judiciously chosen, smaller subset suffices to obtain a set of functions that is sufficiently 
large to separate points of the reduced phase space. 

Let us begin with the holonomies. We will restrict ourselves to three sets of curves: those 
along the R direction of M with oriented length tL q ; those along the equator of S 2 with 
oriented length /x; and those along the longitudes of § 2 also with oriented length /x, where 
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all lengths and orientations are denned using the fiducial triad. (Thus r, /i are positive if 
the tangent to the respective curves are parallel to the triad vectors and negative if they are 
anti-parallel.) Holonomies along these curves suffice to completely determine any invariant 
connection (J2J). 5 

f TL ° TC TC 

h£\A) = exp / dxcr 3 = cos h2r 3 sin— (20) 

Jo 2 2 

hf (A) = exp - ^ dipbn = cos ^ - 2n sin ^ (21) 
Jo 2 2 

hf (A) = exp dfbr 2 = cos ^ + 2r 2 sin ^ . (22) 
Jo z z 

Matrix elements of these holonomies are functions of the reduced connection and constitute 
our configuration variables. Elements of the algebra they generate are almost periodic 
functions of b and c of the form f(b, c) = ^ r f^ T exp | (fib + re), where f^ T G C, /i, r G R 
and the sum extends over a finite set. This algebra is the Kantowski-Sachs analog of the 
algebra of cylindrical functions in the full theory 0, 0, 0]. We will therefore denote it 
by Cylxs- Consider the C* algebra obtained by completing Cylxs using the sup norm. 
The quantum configuration space v4ks is the Gel'fand spectrum of this algebra. From the 
structure of the algebra, it follows that the spectrum is naturally isomorphic to the Bohr 
compactification R| ohr of the Abelian group R 2 J25, 26]. (Recall that in the isotropic case 



the quantum configuration space is isomorphic to Reohr-) The Hilbert space 7i is obtained 
by the Cauchy completion of Cylxs with respect to the natural Haar measure /i on the 
Abelian group R|ohri Bohr' d/i ). Employing the standard bra-ket notation, we can 

define a basis r) in 7i via: 

(b,c\fjt,r) = e^ b+rc) /z,tgR. (23) 

This is an orthonormal basis: 

(fi',T'\fi,r) = 5^^5 T ' T (24) 

where, on the right side, we have the Kronecker symbol, rather than the Dirac delta distri- 
bution. Thus, each basis vector is normalizable and has unit norm. 



In the fully gauge fixed setting, we could omit h$ since it is related to h v by conjugation with exp T3. 
However, we will use a more 'democratic' approach which is applicable also when a is not made to vanish 
by gauge fixing. 
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As one would expect, the configuration — i.e., holonomy — operators hx, hip \ h^ oper- 
ate on 7i by multiplication. However, as in the full theory, these operators are not required 
to be weakly continuous in parameters fi and r, whence there are no operators b, c corre- 
sponding to the connection itself, consequently, although we have only a finite number of 
degrees of freedom, the von-Neumann uniqueness theorem is inapplicable and this quantum 
•;heory is inequivalent to a standard 'Schrodinger quantization' (for a further discussion, see 
26 1). We will see that, as in the isotropic model this inequivalence has important 



25 



consequences. 

For the momentum operators we consider fluxes of triads along preferred 2-surfaces. 
Apart from fixed kinematical factors, they are given by components p&, p c of triads which, 
in view of the symplectic structure (jl2)L are represented by operators 6 

d d 
p b = -z 7 £ pl —, p c =-2ij£ 2 pi —. (25) 



Their eigenstates are the basis states 

Pb\^r) = ^£ Pl fi\fi,r), p c |/i,r) = 74 1 r|/i,r) . (26) 

However, we still have to incorporate the residual gauge freedom which corresponds to 
a parity reflection in the b degree of freedom. Therefore, only those states in TC which are 
invariant under the parity operator fit,: |/i,r) — > \ — fi,r) can belong to the kinematical 
Hilbert space Ti. A basis in TC is thus given by: 

-L[| Ml r> + |-Ai,r>] (27) 

Finally, we express the volume operator in terms of triad operators pb,p c - Recall that the 
region of M under consideration is of the type X x § 2 , where X has length L with respect 
to the fiducial metric. From the classical expression (JT5|) of the volume V of this region, it 
follows that the operator V is given by V = 47r|p b | y^\p c \- It is diagonal in our \fi,r) basis 
and the eigenvalues are: 

V^ = 27T7 3/2 M^M 4i- (28) 
As in the isotropic case, the volume of our cell Ix§ 2 with respect to the fiducial metric can 
be large, but its volume in the 'elementary' state exp | (b + c), is of Planck size. 



6 In this paper, we use the standard quantum gravity convention, £pj = GH. Unfortunately, this is different 
from the convention £p\ — 8irGh used in most of the loop quantum cosmology literature. Hence care 
should be exercised while comparing detailed numerical factors. 
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B. Quantum geometry 



As in the isotropic case the triad operators pb, p c commute whence, in contrast to the 
full theory, one can construct the triad representation. In this sub-section we will study 
the quantum Riemannian geometry using this representation. In the next sub-section we 
will see that the representation is also convenient for dynamics, for it suggests an intuitive 
interpretation for the action of the Hamiltonian constraint. 

Let us then expand a general state in terms of eigenstates of the triad operators, = 
Y^ut Vvl/-*> r )> an d use the coefficients ip^ T to represent the state. This wave function ip^ T 
is supported on the mini-superspace M. coordinatized by pb and p c . In the classical theory, 
already at the kinematical level, the lines pb = and p c = are special: the co-triad (j3J) 
becomes degenerate along pb = (and the line represents the horizon) while it diverges 
on the line p c = (and this line represents the singularity). It is then natural to examine 
wave functions ?/^ T which are supported on lines fi = and r = and ask, already at the 
kinematical level, for the nature of quantum geometry they represent. 

Let us consider the co-triad component uj c = sgn(p c ) \p b \/ ^/\p c \ along dx (see (jlj)). Since 
it involves an inverse power of p c and since there exist normalizable kets r = 0) of the 
operator p c with zero eigenvalue, the naive operator obtained by replacing pb and p c by corre- 
sponding operators is not even densely defined, let alone self-adjoint. To define Cj c , therefore, 
a new strategy is needed. We will follow what is by now a standard procedure, adapted from 
Thiemann's analysis in the full theory j^j]. The first step is to express uj c in terms of the 
elementary variables which do have unambiguous quantum analogs — holonomies and posi- 
tive powers of Pb,p c — and Poisson brackets between them. On the classical phase space, we 
have the exact equality: 

where h x is the holonomy along the interval X, i.e., along the edge of length L (with respect 
to the fiducial metric) we fixed in our construction of the phase space. Then, replacing h x 
and V by their unambiguous quantum analogs and the Poisson bracket by l/ih times the 
commutator, we obtain 

i 



UJr 



2, 7 4, Tr (^fc I - l> 

i ( . c - c c~ . <■ 



( c * c c * c\ 

, f sin — V cos cos -V sin-) . (30) 

V 2 2 2 2 r 



2irj£ P1 
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It turns out that this operator is diagonalized by our basis r) of 

^ C ' M ' ^ = ( y M,r+l - ^T-l) |A*s T ) = — 4l H (V^tTi] ~ ~ 1|) k t) • (31) 

For large r, far away from the classical singularity, these eigenvalues are very close to the 
classical expectation |a> c | = \pb\/ y/\Pc\- This is even true for small /x, i.e. we do not need to 
be far from the horizon. Closer to the classical singularity p c = 0, on the other hand, the 
behavior becomes significantly different from the classical one. At the classical singularity 
itself, where \pb\/ \f\p~c\ diverges, the eigenvalue of u c is zero. 7 

Remark: In the above definition of u c we needed a holonomy in the x-direction. We 
chose to evaluate it using the interval X of length L Q with respect to the fiducial metric. 
While this choice is natural because this interval appears already in the construction of the 
phase space, we could have used an interval of length L S for some 5. Then the resulting 
operator would also have been bounded on the entire Hilbert space H, and |/x, t = 0) would 
again have been eigenvectors of the operator with zero eigenvalue. Thus, the qualitative 
properties of the operator are insensitive to 5. Furthermore, from general considerations one 
can argue P| that 5 should be of the order of 1. However, the value of the upper bound of 
the spectrum does depend on 5. Therefore, as in the isotropic model, the precise numerical 
coefficient in this bound should not be attributed physical significance. 

Let us now turn to the second triad component uj b = sgnp b \/|p c |. Since \f\r~\ is a well- 
defined function of r on the entire spectrum of p c , we can quantize cj& directly. This operator 
is again diagonal in our basis r) and its eigenvalues are given by 

Vb\n, T ) = v^piSS 11 ^) \ZW\ Ia*> t ) ■ ( 32 ) 

Since the operator is quantized directly, the eigenvalues are just the ones one would expect 
from the classical expression (j3J) of u;&. 

The properties of co-triad operators suggest that the lines /x = and r = have very 
different features also in quantum geometry. No quantum effects are manifest at /x = 
which corresponds to the classical horizon. Near and on the line r = which represents the 

7 Exact classical identities such as 1)290 for objects containing inverse powers are available only in homo- 
geneous models. More generally, the bracket between holonomies and volume is more complicated [2^ . 
These complications lead to additional correction terms and the resulting operator is no longer bounded 
on eigenspaces of pb- 
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singularity, on the other hand, quantum effects are large. In particular, they remove the 
classical singular behavior of the co-triad. We emphasize, however, that even in homogeneous 
models, the ultimate test as to whether or not a singularity persists upon quantization can 
only come from studying the dynamics. The key questions are: Is quantum dynamics well- 
defined and deterministic across the classical singularity? And, does this come about without 
significant quantum corrections to geometry near the horizon? Only quantum dynamics will 
tell us if the indication provided by the properties of the co-triad operators are borne out. 

IV. QUANTUM DYNAMICS 

This section is divided into three parts. In the first we present the strategy, in the second 
we obtain the quantum Hamiltonian constraint, and in the third we use it to discuss the 
consequent resolution of the classical singularity. 

A. Strategy 

To bring out the similarities and differences between the reduced model and the full 
theory, it is instructive to begin with the Hamiltonian constraint (|15|) of the full theory: 



We begin by noting some simplifications that arise because of spatial homogeneity. Recall 
that the connection A is related to the spin-connection T defined by the triad and the 
extrinsic curvature K via A = T + ^K. As remarked in section III Al symmetries of the 
model imply that T is the standard magnetic monopole connection; it does not depend 
on the phase space point under consideration. From its expression ©, it follows that its 
curvature Q is given by: 



Since it is a 'c- number' it can be trivially taken over to quantum theory. Now, the curvature 
F of the full connection A can be expanded out as: 



F ab = 2d [a A b] + [A a , A b ] = Q ab + 2 7 <9 [a K 6] + 7 2 [K a , K b ] + 7 [T , K h \ - j[T b , K a \ (35) 




(33) 



Vt = — sin $ T3 d$ A dtp 



(34) 
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Using the expressions (JHJ), © an d Q of T, K and E, it is straightforward to verify that: 

[r, jK] = br 2 cos$d$ Adip = -jdK 

e^K^E-El = (36) 

These equations now imply that the integrand of the Hamiltonian constraint can be written 
simply as 8 

C Ham (x) = e~ l E m E b > (e ijk F* b - 2(1 + i 2 )K\ a K^ = e~ l E ai E b > (e ljk Q k b - 2K [ai K b]j ) . (37) 

Since Q is constant on the entire phase space, the non-trivial part of the Hamiltonian con- 
straint is thus contained just in a term quadratic in extrinsic curvature. 

To pass to the quantum theory, we have to express the right hand side in terms of 
holonomies and triads. There are two possible avenues. 

i) Using (J36)) . the extrinsic curvature terms in the constraint function Cu am (x) can be written 
in terms of curvatures F and Q so that we have: 

C Ham (x) = - 2 e ijk e" 1 E ai E bj [(1 + 7 2 ) Q k ab - F k ah ] . (38) 

The idea now is to use the discussion of co-triads in section IIII Bl to obtain the operator 
analog of e ljk e~ x E ai E h i\ carry to quantum theory trivially and use holonomies to express 
the field strength F as in the full theory. 

ii) Alternatively, since we have completely fixed the gauge freedom to perform internal SU(2) 
rotations, we can regard K itself as a connection. It is a direct analog of the connection A 
in the spatially flat cosmologies where V vanishes. Therefore, we will denote its curvature 
by °F] °F = dK + [K, K]. Then, the constraint functional can also be expressed as: 

C Ham (x) = i e ijh e" 1 E ai E bj [ 7 2 Q k ab - °F k b ] (39) 

Now, one can again use the discussion of co-triads in section IIII Bl to obtain the operator 
analog of e ijk e~ x E ai E b i; carry f2 to quantum theory trivially and use holonomies to express 
the field strength °F. However, now holonomies have to be constructed using the connection 
K rather than A. Since our Hilbert space TC is built from holonomies of A, at first it 

8 This simplification has a natural origin. Because of the right hand side is just tijkE ai E b ^ + F^ b , where 
+ F is the curvature of the self-dual connection + A = T + iK. 
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seems difficult to express holonomies of K as operators on 7i. However, because of spatial 
homogeneity, we can calculate C^amix) at any point. Let us choose a point on the equator and 
use holonomies along the three sets of curves introduced in (J2U|) - (J22|) . Using the expressions 
(J2J) and ( j7| ) of A and K, it follows that, along these curves the holonomies of A and K are 
equal! Therefore, the required holonomies of K can indeed be expressed as operators on TC. 

Both these strategies are viable and, since they begin with exact expressions ()38|) and 
f!39|) of classical constraints, the difference in the resulting operators is just a quantization 
ambiguity. The first strategy is more natural in the sense that it does not require the 
introduction of a new connection. On the other hand, since the new connection K is the 
direct analog of the connection used in spatially flat cosmologies (in particular in [9]), to 
'acilitate comparisons, the second strategy has been used in general, homogeneous models 
sl. l^oj] - For definiteness, we will use the second strategy in the main discussion and comment 
at the end on how the first strategy modifies the final constraint operator. The analysis of 
the resolution of the singularity can be carried out with either methods and the conclusion 
is the same. 



B. Quantum Hamiltonian Constraint 

Starting from the classical Hamiltonian constraint 

C Ham = - 2 jd 3 x e ijh e" 1 E ai E bj [ 7 2 Q k ab - °F k ah ] , (40) 

where the integral is taken over the elementary cell, we wish to pass to the quantum operator. 
The overall strategy is the same as in P, 0, [jjj and the subtleties are discussed in detail in 
jjj. Therefore, here we will present only the main steps and comment on differences from 
the isotropic case treated in [9|. 

In the isotropic case, the elementary cell is a cube of length L with respect to the 
fiducial metric. To express the curvature of the connection and the co-triad operator in 
the Hamiltonian constraint, one uses holonomies along the curves of variable length \i L 
along the edges of this cell. For each value of /i G one obtains a quantum constraint operator 
^Ham- As discussed in |9[, considerations from the full theory imply that jj a should not be 
regarded as a regulator; rather, the /^-dependence of this operator should be regarded as a 
quantization ambiguity. This ambiguity could not be fixed within the reduced model itself. 
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But by making an appeal to results from the full theory, one can fix the value of \i Q using 
the minimum eigenvalue of the area operator in the full theory. (Recall that area enclosed 
by a loop enters while expressing curvature in terms of holonomies.) 

In the present, Kantowski-Sachs model, the overall procedure is the same but minor 
adjustments are necessary because of lack of isotropy. Since S 2 is compact, our elementary 
cell S 2 xl has a geometric edge only along the X direction. As in section IIII Al we will 
supplement it with edges along a longitude and the equator of § 2 . Holonomies along these 
three edges are then given by (|20j) -(|22 JI . Now consider curves of length L Q S along I and of 
length 6 each along the equator and the longitude of S 2 . Then, the co-triad function in the 
Hamiltonian (|3U|) can be expressed as: 

e ljk r l e~ x E^E hk = ^(SttjGC^)- 1 e abc °u k c hf{hf~\ V} (41) 

where on the right hand side the index k runs over x, $, <p and is summed over; hk is the 
holonomy along the edge k; and C x = L Q 5 and = = 5. Note that this is an exact 
equality for any choice of 5. Components of the curvature °F can also be expressed in terms 
of these holonomies. Set 

= hf>hf\hf>)-\ h f)-' (42) 

where k run over x, (p. Then, using (|2~Uj) - (j2~2)) it is straightforward to verify that 



°P**(?)n = " !) + °(( &2 + c2 ) 3/2 ^) ( 43 ) 



where A x $ = S 2 L a = A xip and = 5 2 . For °F X ^ and °F X)tp , is the standard geometric 
relation between holonomies around closed loops, the area they enclose and curvature. For 
°F#,<p, on the other hand, while (J4*H|) continues to hold, the standard geometrical interpreta- 
tion is no longer available because the edges of length 5 along the equator and a longitude fail 
to form a closed loop. Nonetheless, the standard relation is meaningfully extended because 
of spatial homogeneity, S 2 playing the role of the 'effective area' A§^. 

We now have all the ingredients, °F, Q and the triad components, to rewrite the classical 
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constraint ()4Uj) in a way which is suitable for quantization. For the °F-term in (J4*U|) we obtain 
J ^xe i3k °F l ab e~ x E a ^E hk = 2 7 ~ 2 J d 3 x Tr (°F^ n e^e^E^ E bk r l ) 
4(87r 7 3 G)- i y d 3 xdet°^Tr (e^A^hf? - l)£^h { k 5) {h^ 1 ,V}^) + 0((b 2 + c 2 f 2 5) 

+ 0{{b 2 + c 2 f 2 5) . (44) 



Next, using Q = — sin$d$ A dipr^, the Q term in (j4U|) becomes: 

J d 3 xe ijk ni b e~ l E a iE bk = -2 J d 3 xTr (n^n e jM e~ X E al E hk T l ) (45) 

= -8(8tt7G5)- 1 J d 3 xsin^Tr (r 3 C^h^ {hf' 1 , V}) (46) 
= -327r(87r 7 G5)- 1 Tr (r 3 hf{h^- x ,V}) . (47) 

Set 

C^= [-2(7 3 G5 3 )- 1 $>^Tr (hg>^{/i«,y}) ] - [-4(76?^ (r 3 ^{/i?" 1 , V}) ] 

(48) 

Then the classical Hamiltonian constraint is given by 

C Ha m = lim C {S) (49) 

Since all terms in C s are expressed purely in terms of our elementary variables — holonomies 
and triads — which have direct operator analogs, passage to quantum theory is now straight- 
forward. We obtain: 

\ ijk 

w s-i fn - fib 5b . 5c 5c ( . 5b~ 5b 5b~ . 5b s 

= 4? 7 L, ) 8 sin — cos — sin — cos — sin — V cos cos — V sin — 

w PV \ 2 2 2 2 \ 2 2 2 2 / 

( . . 2 5b 2 5b 2r2 \ / . 5c- 5c 5c- . 5c\\ 
+ I 4 sin — cos h 7 I I sin — 1/ cos cos — V sin — 1 1 (50; 

v 2 2 y v 2 2 2 2 y y 

(This is a special case of Eq. (26) in j30|.) 

The action of this operator on the eigenstates r) of pb and p c is given by 

#V r) = (2 7 3 5 3 4 1 )" 1 [ 2 ( W " U) (51) 
x + 25, r + 25) - |/i + 25, r - 25) - |/x - 25, r + 25) + |/i - 25, r - 25)) 
+ (V; t , T+5 - ^-,0(1// + 45, r) - 2(1 + 2 7 2 5 2 ) r) + |/x - 45, r))] . 

20 



However, while the classical constraint is a real function on the kinematical phase space 
r, C 5 fails to be self-adjoint on the kinematical Hilbert space 7i. We therefore need to add 
to it its Hermitian adjoint. 9 The result is an operator I := \{C^ + C^) given by: 

C£L\V,T) = (2 7 3 5 3 4l) _1 [{V^ + s >T -V^S ) r + V^ +3 s,r+25-V^ +S)T+2S )\fX + 25,T + 25) 
-{y^+5,r - V»-S,t + ^t+3<5,r-25 ~ ^t+(5,r-2(j) |# + 25, T - 25) 
— {Vh+5,t - y^-5,r + V^-S, T +2S - V^-35,r+2<y) \fi> - 25, t + 25) 
+ (V M+ 5,r - Vfi-8,r + V^-S, T -2S - X^u-35,-r-2<y)|A* 25, T 25) 

+|C^,t+5 - V^ )T _5 + Vp+^T+s - Vfj, + 4s,T-s)\fJ< + 45, r) 
-(1 + 2 7 2 5 2 )(^, r+(5 - V^ s )\fi, t) 

+ U V v,r+s ~ y^,r-8 + V^4S,r+s ~ V^_45, T -*)|^ - 45, r)] . (52) 

Physical states in quantum theory are those which are symmetric under the 'parity op- 
erator' fib and lie in the kernel of the operator Cgfl v . That is, the only non-trivial quantum 
Einstein's equations are: 

(vi>|n b = o , (*|c£ v = o (53) 

where is an element of the dual Cyl^ s of the space Cylxs of finite linear combinations 
of almost periodic functions of b, c. Let us expand (\1/| using eigenbras (//, r| (in Cyl^ s ) of 
the triad operators: 

(*| = X>*r(/i,r| (54) 

fJ,,T 

Then, ij;^ can be regarded as wave functions in the triad (or, Riemannian geometry) repre- 
sentation. To exhibit the action of Cgfl v on t/> M)T we are led to separate the cases /x > 45 from 
ji < 45 because of the absolute- value in the volume eigenvalues. For \i > 45, quantum 



9 Although the self-adjoint form of the Hamiltonian constraint was discussed briefly in [^3, this point was 
ignored in the detailed treatment of M. A detailed treatment of the self-adjoint constraint and its semi- 
classical implications can be found in fed . It is used in a crucial way to obtain the physical Hilbert 
space in lid). 
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Einstein's equation becomes 

-2<5,r+2<5j 

+(\/|t + <5| - y/\r-5\) ((/i + 25)^+4<5,r - (1 + 27 2 «5 2 )/x^,r + - 2<J)V M -«,r) 

+25(a/|t - 2<5| + a/|t[) (V> m -2£,t-2<5 - l(>n+28,T-26) 

= 0. (55) 
For ji — 0, we have 

V|r + 2<J|V^,r+a« + + - V|r - <5|)^,r - y/fr ~ ^\ihe, T -26 = ; (56) 

for n = 5: 

2(v/|r + 25| + v^i)fe,r + 25 + 2(^1^ + 251 - vW)Vv+2<5 
+ + - y/\r-5\) (3fe, T + 2^ 35l r - (1 + 2 7 V)Vv) 
-2(v/|r-25| + x/M)^,^ - 2(^-251 - vH)Vm = > ( 57 ) 

for /i = 25: 

(a/|t + 25| + a/M)^45,t+25 - a/|t|^o, t+ 2 ( 5 

+ (^ + 51 " Vk-^l) ( 2 ^,r + - (1 + 2 7 2 5 2 )^,r) 

- ( y/\T-25\ + V^)^45,r-25 + W\*Po,T-25 = , (58) 

and for fi = 35: 

2(V|r + 25| + \/jrf)(^5*,r+25 - ^<5, 

+( V / k + ^l " (5^75,r + 2?/V - 3(1 + 2-f8 2 )ip 3S ,r) 

-2(y/\T - 25\ + V^)(fe,r-25 - At-25) = . (59) 

For values of /i not an integer multiple of 5 one can derive similar expressions, but they 
will not be used in what follows. As remarked above, strictly speaking, the role of (J55J) 
is only to select physical wave functions. However, intuitively it can also be thought of 
as the quantum evolution equation. Recall that p c can be regarded as an internal time in 
the classical theory of the model. Therefore, in the quantum theory, one may regard r as 
internal time. Then, ()55|) can be interpreted as the quantum Einstein's equation, which 
evolves the wave functions in discrete steps of magnitude 25 along r. 
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So far, the parameter 5 — and hence the size of the 'time-step' — is arbitrary. In the 
classical theory, we recover the Hamiltonian constraint only in the limit 5 goes to zero. This 
is because, while the expression (j3J) is exact for all loops, curvature °F can be expressed 
in terms of holonomies only in the limit in which the areas of loops shrink to zero (see 
f|43)l ). In quantum theory, however, the straightforward limit lim^o diverges (for the 
same reasons as in the isotropic case Q]). The viewpoint is that this is occurs because the 
limit ignores the quantum nature of geometry, i.e., the fact that in full quantum geometry, 
the area operator has a minimum non-zero eigenvalue. The 'correct' quantization of the 
Hamiltonian constraint Cnam nas to take in to account the quantum nature of geometry. 

To do so, note first that the holonomies used in the expression of °F define quantum 
states with fi = 5 and r = 5. Using these states, we can calculate the quantum area of each 
of the three faces of the elementary cell § 2 x X, enclosed by the three curves — one along 
x, one along a longitude and one along the equator of § 2 . The area operators defined by 
faces S x ^,S X)ip and S^ )ip are, respectively, 27r|p&|, 2tt | j3t.| and 7r|j3 c | (see fjl4|) ). The quantum 
geometry states defined by the three holonomies are eigenstates of these area operators. 
Furthermore, they have the same eigenvalue: wySi^. Now, in the full theory, we know that 
the area operator has a minimum non-zero eigenvalue, a a = 2\^3irj£pi. The viewpoint, as 
in P|, is that in the calculation of the field strength °F, it is physically inappropriate to try 
to use surfaces of arbitrarily small areas. The best one can do is to shrink the area of the 
loop till it attains this quantum minimum a Q . This implies that we should set 5 = 2\/3. To 
summarize, by using input from the full theory, we conclude that the quantum Hamiltonian 
constraint is given by C {So) with 5 Q = 2^3. 

Remarks : 

i) Note that the character of the difference equation changes depending on whether fi > 45 
or fi < 45. In particular, for \x > 45, the knowledge of the wave function at 'times' r + 25 
and r determines only the combination 1^^-25^-25 — Vv+2S,r-2<5 at 'time' r — 25 via (J33j). On 
the other hand, for \x = the wave function at r + 25 and r determines ip26,r+2S completely 
via (JHEJ). We will return to these differences in section ITVCl 



ii) As discussed in 



19, 
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34j in detail, one can recover the Wheeler-DeWitt equation 



with a specific factor ordering by taking a systematic limit of quantum Einstein's equation. 
In the present model, the difference equation reduces to the following differential equation 
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for the wave function ^(pb,p c ) = ty'ivbhP ,vch£ % x 011 the classical minisuperspace M.: 

pi yf c dpbdpc 4y/pZ dpi 2^/p~ c dp b ) yfp c } 

For /i 5, r ^> 5, solutions of the discrete equation (J55j) can be approximated by those of 
this Wheeler DeWitt equation in a precise sense. In essence, this Wheeler-DeWitt equation 
is obtained from by ignoring quantum geometry effects, i.e., by an appropriate 5 — > 
limit. For the approximation by a differential equation, higher derivatives of the wave 
function ip must be sufficiently small. Thus, the limit 5 — ► for the equation (}55|) exists 
only under additional assumptions on the wave function. It does not exist for the operator 
C'grL by itself. 

C. Absence of Singularity 

Recall that the classical singularity occurs at p c = 0. The Wheeler DeWitt equation (JfiUJ) 
is manifestly singular there. One can multiply it by ^Jp~ c and define a new 'internal time, 
p c = \np c and obtain a regular equation in variables Pb,p c - However, since p c — > —00 at 
the singularity, this regular equation does not let us evolve across the singularity. Thus, the 
overall situation is the same as in the isotropic model. 

To obtain the Wheeler-DeWitt equation from the 'fundamental' difference equation (|55j) . 
we had to ignore quantum geometry effects. Thus, the key question is whether this failure 
of the Wheeler-DeWitt equation is an artifact of the approximation used. Can quantum 
geometry effects make a qualitative difference as they do in the isotropic model? We will 
now argue that the answer is in the affirmative. 

Let us analyze the 'evolution' given bv the difference equation ()55jl . Does this evolution 
stop at t = 0? As in the isotropic case P,E3, we wm use ^ ne quantum constraint (|HH)l as 
a recurrence relation starting at positive r and evolve toward smaller values. However, now 
a new twist arises because for generic values fi > AS, determines only the difference 
^h+25,t-25 ~ i>n-2S,T-26 as a function of the initial values of the wave function. Therefore, 
quantum Einstein's equation has to be supplemented with an appropriate boundary condi- 
tion. From (jf)5j) which holds for fi > 45, one might conclude that, even if we restrict ourselves 
to the 'lattice' fi = n8, one would have to specify the wave function at /x = 0, 5, 25, 35 at 
each 'time step'. However, as remarked at the end of section HVBl the form of the difference 
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equations for // < 45 is different. An examination of ()56j) and (joT|) reveals that it is in fact 
sufficient to specify the wave function ip^ T just at fi = and = 5 at each 'time-step'. This 
is a mathematically viable choice of the boundary condition and heuristically it corresponds 
to providing data the 'horizon'. It will turn out that the issue of the resolution of singularity 
is insensitive to the precise choice of the boundary condition. We will therefore postpone 
the discussion of physically appropriate choices until after the main discussion. 

Let us then start at some finite positive r and carry out a backward evolution using (J57)|) . 
In contrast to the Wheeler-DeWitt equation ()60|). the coefficients in the difference equation 
()55j) are always regular. At first one might think that this regularity by itself would be 
sufficient to ensure that the evolution would be well-defined across the singularity. Note 
however that this need not be the case. For, the backward evolution continues only as long 
as the coefficient of the wave function at r — 25 is non-zero. The key question therefore is 
whether this coefficient vanishes. The issue is subtle. For instance, the coefficients in the 
non-self-adjoint operator (C s y are also non-singular. However, if one uses it in place of 
C^ am , one finds that the coefficient vanishes and one can not evolve across the singularity. 

For our quantum Einstein's equation (jHHjl . this coefficient is given by y/\r — 25\ + \f\r~\- 
By inspection, it never vanishes. Therefore, the backward quantum evolution remains well- 
defined and determines the wave function not only for r > but also in the new region 
with r < 0. In this precise sense, the classical black hole singularity can be traversed using 
quantum evolution and thus ceases to be a boundary of space-time. 

Next, recall from section II V Al that we could have begun with the expression (J38j) of 
the Hamiltonian constraint in the classical theory and then proceeded with quantization. 
What would be the status of the singularity resolution with this choice? The procedure 
to construct the quantum Hamiltonian constraint would be identical. However, the final 
result would be slightly different: 7 2 in (|B3j) would be replaced by 1 + 7 2 . That is, in the 
new quantum Einstein equation, the coefficient of ip^, r in (|B3|) would be modified and the 
equation would become: 



+(V\r + 5\ ~ V\t-S\) ((/! + 25)^ +4 5,r - (1 + 2( 7 2 + l)5 2 )/i</v + - 25)^_ 4(S , r ) 




V-2<5,t-2<5 — '4 ) tM+2S,T-2S) 







(61) 
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(The same replacement holds for /i < 45.) Since the coefficient of VV-2<5,t-2<5 — Vv+2<5,t-2<5 is 
unaltered, one would again be able to 'evolve' across the singularity. Thus, the conclusion 
is robust within quantum geometry. 10 

To conclude, we will briefly return to the question of the boundary condition that are 
needed to make the evolution well-defined even away from the singularity. As noted earlier, 
in the backward evolution given by ()56|) . ip28,r is determined by values of the wave function 
at t > t , so long as for r D 7^ 0. Furthermore, ()55|) then determines the wave function at 
U = (2 + 4n)5 for all n G N and at that r . One could now use pre-classicality arguments 
|35j ] (i.e. require that the wave function not oscillate on small scales) at large /i to choose the 
boundary values, eliminating entirely the need for specifying boundary conditions. While 
this could be reasonable in semiclassical regimes, for small |r| it would be questionable 
to use pre-classicality in the /i-direction even if a pre-classical solution exist (which is not 
guaranteed, see e.g. the analysis in 3, l^ll)- More importantly, the slice r = is special 
because ip2S,o is not determined through ()56|) or otherwise. At this point, the condition (JoT)j) 
evaluated at r = —25 gives instead a condition for previous values of the wave function, 
translating into one condition for the initial values. At r = 0, however, ip28,o has to be 
specified as boundary value in addition to V'o.o because, unlike similar situations in isotropic 
models, it does not drop out of the evolution: it is needed in the recurrence (jo^j) for n = 45, 
t = 25. There is thus no reduction in conditions on the wave function from the fact that 1^26,0 
drops out of (JoT)|) ; a condition is simply transferred from boundary values to initial values. 
Note however that, even though the line r = (corresponding to the classical singularity) 
does show special behavior, evolution does not break down there; it is only the boundary 
value problem which changes character. 

While such a boundary value problem is mathematically well-defined, the resulting theory 
is not necessarily physically correct. For, in a physically interesting theory one would expect 
that, well away from the singularity, there should exist semi-classical solutions which are 
peaked on the classical trajectories. From Fig. 1 it is clear that the semi-classical solution 
peaked on the trajectory labelled by mass m would be sharply peaked at r = ±4m 2 on 
the line /1 = 0, whence semi-classical states peaked at different classical solutions will have 



10 However, the relation between the wave function and the Wheeler-DeWitt approximation is less direct 
with this form of the constraint. 
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quite different forms near fi = 0. Therefore, the theory obtained by simply fixing the wave 
function on (or near) jj, — is not likely to admit a sufficiently rich semi-classical sector. 
Indeed, it is not clear that any boundary condition at (or near) /i = will be physically 
viable. Rather, one may have to impose the boundary condition at fi = oo, e.g., by requiring 
that the wave function should vanish there. Indeed, the form of dynamical trajectories of 
Fig. 1 implies that every semi-classical state would share this property. Moreover, the 
corresponding condition does hold in the closed isotropic model with a massless scalar field 
as source 3sj|. As in that case, one would expect that these semi-classical states would span 
a dense subspace in the physical Hilbert space. If so, requiring that the wave functions 
vanish at ji = oo would be physically justified. Heuristics suggest that this strategy is viable 
for the Wheeler-DeWitt equation (J6U|) . However, detailed numerical analysis is necessary to 
establish that the difference equations (jo3j) or (j61|) admits solutions satisfying this condition 
(or a suitable modification thereof) for all r and that its imposition makes the evolution 
unambiguous. 



V. DISCUSSION 



Results of the last two sections support a general scenario that has emerged from the 
analysis of singularities in quantum cosmology. It suggests that the classical singularity 
does not represent a final frontier; the physical space-time does not end there. In the Planck 
regime, quantum fluctuations do indeed become so strong that the classical description 
breaks down. The space-time continuum of classical general relativity is replaced by discrete 
quantum geometry which remains regular during the transition through what was a classical 
singularity. Certain similarities between the Kantowski-Sachs model analyzed here and a 
cosmological model which has been studied in detail [l^ suggest that there would be a 
quantum bounce to another large classical region. If this is borne out by detailed numerical 
calculations, one would conclude that quantum geometry in the Planck regime serves as a 
bridge between two large classical regions. Space-time may be much larger than general 
relativity has had us believe. 

However, as indicated at the end of section HV Bl significant numerical work is still needed 
before one can be certain that this scenario is really borne out in the model. Moreover, this 
is a highly simplified model. It is important to check if the qualitative conclusions remain 
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robust as more and more realistic features are introduced. First, one should extend the 
analysis so that the space-time region outside the horizon is also covered. A second and 
much more important challenge is incorporation of an infinite number of degrees of freedom 
by coupling the model, e.g., to a spherically symmetric scalar field. First steps along these 

|40| and one can see that the evolution still extends beyond the 



lines have been taken 



classical singularity |4l|. But a comprehensive treatment still remains a distant goal. 
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